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Abstract 

We investigate the connection between local and global dynamics in 
the Fermi - Pasta - Ulam (FPU) {3 - model from the point of view of 
stability of its simplest periodic orbits (SPOs). In particular, we show 
that there is a relatively high q mode (g = 2(A'^ + l)/3) of the linear lat- 
tice, having one particle fixed every two oppositely moving ones (called 
SP02 here), which can be exactly continued to the nonlinear case for 
Af = 5 + 3m, m = 0, 1, 2, . . . and whose first destabilization, E2u, as the 
energy (or /3) increases for any fixed A^, practically coincides with the onset 
of a "weak" form of chaos preceding the break down of FPU recurrences, 
as predicted recently in a similar study of the continuation of a very low 
(g = 3) mode of the corresponding linear chain. This energy threshold per 
particle behaves like oc A'^"^. We also follow exactly the properties 
of another SPO (with q = (A^^ + l)/2) in which fixed and moving parti- 
cles are interchanged (called SPOl here) and which destabilizes at higher 
energies than SP02, since oc A'^"^. We find that, immediately after 
their first destabilization, these SPOs have different (positive) Lyapunov 
spectra in their vicinity. However, as the energy increases further (at fixed 
A'^), these spectra converge to the same exponentially decreasing function, 
thus providing strong evidence that the chaotic regions around SPOl and 
SP02 have "merged" and large scale chaos has spread throughout the 
lattice. Since these results hold for A'' arbitrarily large, they suggest a 
direct approach by which one can use local stability analysis of SPOs to 
estimate the energy threshold at which a transition to ergodicity occurs 
and thermodynamic properties such as Kolmogorov - Sinai entropies per 
particle can be computed for similar one - dimensional lattices. 



1 Introduction 

The transition to widespread chaos in Hamiltonian systems of many degrees 
of freedom has been the subject of intense investigation for more than fifty 
years, see e.g. 0; |2 S 01- It received great impetus following the pioneering 
work of Fermi, Pasta and Ulam 0, who were the first to study thermalization 
in one - dimensional lattices of N particles, with linear and nonlinear nearest 
neighbor forces, as a parameter multiplying the nonlinear terms in the equations 
of motion becomes greater than zero. Surprisingly, they observed that when this 
parameter is relatively small, energy cquipartition docs not occur even after 
very long integration times, as only a small number of (low q) modes of the 
corresponding linear lattice recurrently exchange the total energy among them. 
Of course, when the nonlincarity or the energy exceeds a certain threshold, these 
so called FPU recurrences break down, large scale chaos prevails and a type of 
ergodicity sets in, whereby almost every orbit explores almost all of the available 
phase space of the system. 

One of the first attempts to explain this phenomenon is due to Izrailev and 
Chirikov Q , who argued that the breakdown of FPU recurrences is related to 
the overlap of major resonances known to lead to large scale chaos in N - degree 
of freedom Hamiltonian systems 0. As it was later realized, however, a weaker 
form of chaos caused by the interaction of the first few FPU modes (with low 
q in Fourier space) appears to be sufficient for cquipartition among all modes 
to occur 01 . Finally, very recently, Flach and coworkers ^3 discovered that 
this transition to the so called "weak" chaos, in fact, coincides with the first 
destabilization of one of the lowest {q = 3) normal mode of the linear lattice, as 
it is continued by increasing the nonlincarity parameter. In fact, their results 
apply more generally to the lowest q <ti N modes which also turn out to be 
highly localized in q - space. These nonlinear modes represent examples of 
what we call simple periodic orbits (SPOs), where all particles return to their 
initial condition after only one oscillation, i.e. all their mutual rotation numbers 
are unity [llj. Such lowest q mode SPOs, were actually termed q - breathers 
due to their exponential localization in Fourier space |lfll |. 

In this paper, we investigate further the connection between local and global 
dynamics of the FPU lattice, by studying the stability properties of its SPOs. 
In particular, we consider the FPU Hamiltonian 



often called the FPU - model, as it only contains the term with quartic 
nonlinearities. The Xj represents the displacement of the jth particle from 
its equilibrium position, Xj is the corresponding velocity, /3 is a positive real 
constant and E is the total energy. As with the original FPU problem, we will 
concern ourselves only with fixed boundary conditions, whereby particles with 
index j = 0, iV + 1 arc stationary for all time. 

In particular, we will examine an SPO which keeps every third particle fixed. 
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while the two in between are performing exactly opposite motions. This mode 
was studied in where its stability was analyzed by means of Mathieu equa- 
tions and has also been discussed in in connection with the occurrence of 
sinusoidal waves in nonlinear lattices. This solution, called SP02 from here on, 
will be compared to an orbit we call SPOl, which keeps every two particles fixed, 
with the ones in between executing exactly opposite oscillations. This latter one 
was originally mentioned in a paper b y O oyama et. al. |l4l| and later studied 
analytically by Budinsky and Bountis [l5|| to determine its stability properties 
in the thermodynamic limit of large TV and E with ^ fixed. Recently, it was 
revisited by Antonopoulos et al. in a study of diiferent SPOs and different 
Hamiltonians, from the viewpoint of connecting their local and global stability 
properties. 

Our first result about the SP02 orbit is that the energy per particle of its first 
destabilization goes to zero faster than SPOl, by a law oc in contrast 
to the SPOl orbit whose law is cx as iV -> oo lll[ll|. This impKes 
that if chaos is to spread in the nonlinear lattice as a result of the destabilization 
of SPOs, it might be more useful to look closely at the properties of SP02, as 
that becomes unstable much earlier than SPOl, as N increases. 

Remarkably enough, when we do this we discover that the energy (or (3) val- 
ues of the first destabilization as a function of N practically coincide with those 
found by 13 01 for the transition to "weak" chaos and for the destabilization 
of the q = 3 mode. Our numerical results and their analytical formula are in 
excellent agreement. 

We then examine the dynamics in more detail, following the first destabiliza- 
tion of SPOl and SP02 at E = Eiu and E = £'211 respectively, at high enough 
N . In particular, we choose initial conditions in the vicinity of these orbits and 
find that the (positive) Lyapunov exponents, at energies just above Eiu and 
i?2u) fall off to zero following distinct curves, both for SPOl, which destabilizes 
by a period - doubling type of bifurcation and SP02, which exhibits a complex 
instability with its monodromy eigenvalues exiting the unit circle in complex 
conjugate pairs. However, as the energy increases further, the Lyapunov spec- 
tra near SPOl and SP02 begin to converge, at some E > Eiu > £^2ti, to the 
same functional form, implying that the chaotic regions of SPOl and SP02 have 
"merged" and large scale chaos has spread in phase space. 

The function to which the Lyapunov spectra converge is a nearly exponen- 
tially decaying curve of the form 

L,{N)o^e-"^,i = l,2,...,K{N) (2) 

at least up to K{N) sa as we have discussed at length in a recent publication 
[lll |. This function provides, in fact, an invariant of the dynamics, in the sense 
that, in the thermodynamic limit, we can use it to evaluate the average of the 
positive Lyapunov exponents (i.e. the Kolmogorov - Sinai entropy per particle) 
and find that it is a constant characterized by the value of the exponent a 
appearing in 

Thus, we argue that studying the local dynamics around some of the sim- 
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plest periodic orbits which destabilize at low energies, opens a "window" into 
the "global" dynamics of nonlinear lattices. Furthermore, by computing and 
comparing Lyapunov spectra in their vicinity, it is possible to gain valuable in- 
sight into the conditions for large scale chaos and ergodicity, so that we may be 
able to define probability distributions and compute thermodynamic properties 
of the lattice, as E and N increase indefinitely with E/N fixed. 

Our paper is organized as follows: In section [21 we provide analytical expres- 
sions of the SPOl and SP02 solutions under study and describe in detail their 
stability properties for an arbitrarily large number of particles of the FPU /3 - 
model with fixed boundary conditions. Comparing with similar findings in the 
literature, we observe that our results accurately predict the onset of "weak" 
chaos preceding the break down of FPU recurrences, although the reasons for 
this agreement are still under investigation. In section [SJ we use our results 
on the convergence of Lyapunov spectra, as the energy is increased beyond the 
destabilization thresholds of SPOl and SP02 to estimate the onset of large scale 
chaos and thermodynamic behavior in the lattice and in section^lwc present our 
conclusions. Wc thus believe that today, 50 years after its famous discovery, the 
Fermi - Pasta - Ulam problem and its transition from recurrences to globally 
chaotic behavior is still very much alive as a topic of intense research into some 
truly fundamental questions connecting classical and statistical mechanics [l^ . 



2 Simple Periodic Orbits and Stability Analysis 
2.1 Analytical results for SPOl 

Let us start by describing briefly some analytical results concerning SPOl, as 
this particular mode has been studied recently by Antonopoulos et al. [lH and 
also previously in 0; 0| . 

We consider, for this reason, a one - dimensional lattice of N particles with 
equal masses and nearest neighbor interactions with quartic nonlincaritics (/? 
- model) which is given by the FPU Hamiltonian 0}, with fixed boundary 
conditions 

xoit) =XN+i{t) = Oyt- (3) 

For P — 0, Hamiltonian describes a system of coupled harmonic oscilla- 
tors and hence all solutions can be written as combinations of independent 
normal modes whose individual energies are constant in time. Since, in that 
case, the spectrum has frequencies ujq = 2sin(7r(j'/2(A^-t- 1)) which are rationally 
independent, all solutions are quasiperiodic, and hence the only strictly per iodic 
solutions are the normal modes, with frequencies ujq, q = 1,2, . . . , N [allOtllq. 
That these modes can be continued for /3 > is a consequence of a famous 
theorem by Lyapunov , based on the assumption that no ratio of linear fre- 
quencies LOq/uJr is an integer, for q,r = 1,2, . . . , N , which holds in this case. 
These solutions are examples of what we call Simple Periodic Orbits (SPOs), in 
which all particles return to their starting point after one maximum (and one 
mininmm) in their oscillation \ll\ . 
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Let us consider one such orbit - we shall call SPOl - which is specified by 
the conditions 



iV - 1 

X2j{t) ^ 0, X2j-l{t) = -X2j + l{t) = x{t), i = 1, . . . , (4) 

and exists for all odd N, keeping every even particle stationary at all times. It 
is not difficult to show that this is, in fact, the q = {N + l)/2 mode of the linear 
lattice with frequency ujq = \f2. The remarkable property of this solution is that 
it is continued in precisely the same spatial configuration in the nonlinear lattice 
as well, due to the form of the equations of motion associated with Hamiltonian 

o, 

(0 = - 2xj + Xj^i +f3(^{xj+i ~Xj)^ - {xj -Xj^i)^^, j = 1, . . . ,iV (5) 

which reduce, upon using with (PJ to a single second order nonlinear differ- 
ential equation for x{t), 

i(t) = -2x{t) - 2/3x^(t) (6) 

describing the oscillations of all moving particles of SPOl, with j = 1, 3, 5, . . . , A^. 
For the stationary particles j = 2, 4, 6, . . . , TV— 1 of course, we have x{t) = 0, Vt > 
0. The solution of @ is well known in terms of Jacobi elliptic functions ^3 
and can be written as 

x{t) =Ccii{Xt,K^) (7) 

where 

^ ^ (3{1-2k')' ^ ^ l^2n^ 

and is the modulus of the en elliptic function. The energy per particle of 
SPOl is then found to be 

^ - lc'(2 + C'3) - '''^^^""^^^ (Q) 

by substituting simply the solution x{t) of (0 in Hamiltonian Q). 

The linear stability analysis of the SPOl mode is straightforward and was 
carried out recently in [lH using Lame equations, Hill's determinants and Flo- 
quet theory. Plotting the first destabilization energy for this orbit as a function 
of N with solid lines in Fig. ^ we observe that the corresponding energy density 
threshold decreases with N following a simple power law (x 1/N (dashed 
line). Following such an approach, we find, for example, for (3=1 and = 11, 
that SPOl destabilizes for the first time when w 1.93. 

2.2 The solution of SP02 

Let us now turn to the second simple periodic orbit studied in this paper, which 
we call SP02. We impose again fixed boundary conditions and consider lattices 
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Figure 1 : The solid curve is the energy per particle of the first destabilization 
of the SPOl for /3 = 1 while the dashed line is the function oc jj-. 



consisting of = 5 + 3m, m = 0, 1, 2, . . . particles, in which every third particle 
is fixed, while the two in between move in opposite directions as follows 

iV - 2 

X3jit) = 0, j = f,2,3...,^— , (10) 

Xj{t) = ~xj+i{t)^x{t), j ^l,A,7,...,N-l. (11) 

This solution, in fact, corresponds to the q = 2{N + l)/3 normal mode with 
frequency ujq = a/3 of the linear system and can also be continued in exactly the 
same form in the nonlinear case /? > 0, due to the symmetry of the equations 
of motion, 

Xj{t) ^ x^+i-2xj+Xj^i+(3(^{x.i+i~Xjf-{xj~Xj^ify j ^1,...,N (12) 

which, under the above conditions, H1U|) . collapse to a single second order 
nonlinear differential equation very similar to (jSj) 

i(t) = -3x{t) -9f3x^{t). (13) 

As before, this equation describes the moving particles of the lattice, while 
the stationary ones satisfy x{t) = 0,Vt > 0, for j = 3, 6, 9, . . . , iV — 2. The 
solution of equation H13|) is again given in terms of the Jacobi elliptic functions 
[l^ and is written in the form 

x{t) =Ccii{Xt,K^) (14) 

where 

= 3/5(1-2^2)' ^ (15) 
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and is, again, the modulus of the cn elliptic function. The energy per particle 
of the SP02 mode is found to be now 

E _ 2k^{1~k^) 
N + 1 ^ 3/3(1 - 2k2)2 

by simply substituting the solution x{t) of ifT^ in Hamiltonian QJ. 

In order to perform the linear stability analysis of the SP02 mode wc set 
Xj = Xj + Uj in the equations of motion H12(l and keep up to linear terms in the 
small displacement variable yj. We thus get the variational equations for this 
orbit in the form 

yj{t) = Agj/j-i + Aivj + A2yj+u j = 1, 4, 7, . . . , iV - 1, (17) 
yj{t) = A22/j„i+Aiy, +A3%+i, J = 2,5,8,...,7V, (18) 
yj{t) = A3(j/,-i-22/,+yj+i), J =3,6,9,...,7V-2 (19) 

where y^ — yN+i — and 

Ai = -2- 15/3*2 (i), (20) 

A2 = l + 12Px^{t), (21) 

^3 = l + 3Px^{t). (22) 

Unfortunately, it is not as easy to uncouple the above linear system of 
differential equations and study the stability of the SP02, in terms of inde- 
pendent Lame equations, as we were able to do with SPOl We can, 
however, compute numerically with arbitrary accuracy and for every given 
N ~ 5 + 3m, jTi = 0, 1, 2, . . . the complex eigenvalues Ai, z = 1, . . . , 2N of the 
corresponding monodromy matrix and characterize the stability of the SP02 by 
their position on the complex plane with regard to the unit circle. 

We have thus computed, for many values of = 5 + 3m, m = 0, 1, 2, ... , 98, 
the energy E2u{N) of the first dcstabilization of the SP02 for /3 = 1 and have 
plotted the results with solid lines in Fig. |2Ia). As wc sec, the energy density 
at the first instability decreases following a power - law cx (dashed 
line) which is faster than the SPOl solution wc discussed earlier, see Fig. ^ 
Following this approach, wc find, for (3 = 1 and = 11, that SP02 destabilizes 
for the first time when E2U ~ 0.153. 

Interestingly enough, if we calculate the eigenvalues of the monodromy ma- 
trix of the SP02 for greater energies, we find that it becomes again stable, 
beyond a new critical energy E2s{N). In Fig. El^b) we plot this restabilization 
energy density of SP02 as a function of log N and observe that it approaches a 
constant as N tends to infinity. 

Pursuing further this restabilization phenomenon, wc have estimated the 
"size" of the islands of stability around SP02 for energies above the energy 
threshold E^MV), using the method of the Smaller ALignmcnt Index (SALI) 
[19I : Eot I211 12^ . This index has proved to be very efficient for distinguishing 
rapidly and with certainty regular vs. chaotic orbits, as it exhibits completely 
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Figure 2: (a) The solid curve is the energy density of the first destabihzation 
of the SP02 obtained by the numerical evaluation of the eigenvalues of the 
monodromy matrix, while the dashed line is oc (b) The energy density 
of the first restabilization of the SP02 obtained by the same method as in (a). 
In this figure (3 = 1. 



different behavior in these two cases: It fluctuates around non - zero values 
for regular orbits, while it converges exponentially to zero for chaotic orbits. 
SALI is particularly useful in the case of many degrees of freedom, where very 
few such methods are available beyond the cumbersome and often inconclusive 
calculation of the maximal Lyapunov exponent. 

We thus observe the following: As the energy E increases beyond the resta- 
bilization threshold E2s{N), for fixed A^. the "size" of the island around SP02 
changes very little, compared with the growth of the system's available phase 
space. Moreover, if wc keep ^ fixed, thus holding the "radius" of the en- 
ergy surface nearly constant, we find that the "radius" of the SP02 island 
diminishes as a function of E. This is done by changing one of the par- 
ticles' position and momentum by Ax and Apx away from its SP02 values 
while keeping the energy constant and using SALI to estimate Axmax which 
we reach chaos (e.g. with ^ = 4 and A^ = 5,8,11,14, we find respectively 
Axmax ~ 0.01,0.005,0.0025,0.002). We, therefore, conclude that the island 
of ordered motion around the SP02 solution should be of no consequence to 
the statistical properties of the lattice, such as crgodicity and the definition of 
thermodynamic quantities. 

Moreover, as we observe in Fig. Efa), (c), the kinds of bifurcation leading 
to instability for these SPOs are very different: In the case of SPOl, Fig. Eta) 
shows that the bifurcation is of the period - doubling type, as one pair of real 
eigenvalues is seen to exit the unit circle at -1, while Fig. |3fc) shows that we have 
complex instability in the case of SP02. This is also indicated by the positive 
Lyapunov exponents in the immediate vicinity of the SPOs (about 10^^^ from 
them), which are closely connected with the eigenvalues of the corresponding 
variational equations and are shown in Fig. Ol^b), (d) with solid lines. Of 
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course, moving away from the two modes (within a range of about 10~^^ - 
10~^), the Lyapunov spectra change into the famihar form of two smoothly 
decaying, evidently different curves, plotted with dashed hues in Fig. Elb), (d). 
This fact suggests that the chaotic regions near these modes arc separated from 
each other in phase space. 

It is quite interesting to observe that the chaotic region about an unstable 
SPO can be isolated in phase space from the chaotic motion occurring in differ- 
ent domains. In fact, there may be several such domains embedded into each 
other. For example, in the = 5 FPU (3 - model, when SPOl becomes un- 
stable, a "figure eight" chaotic region becomes clearly visible in its immediate 
neighborhood, on a Poincare surface of section (xi,ii) taken at times when 
^3 = (see Fig. 0}. Even though the SPOl mode is unstable, nearby orbits 
oscillate about it for very long times, forming eventually the "figure eight" we 
see in the picture. 

More surprisingly however, starting at points a little further away, a different 
chaotic domain is observed which bears a vague resemblance to the "figure eight" 
and does not spread to the full energy surface. Of course, if one chooses more 
distant initial conditions a large scale chaotic region becomes evident on the 
surface of section of Fig. 0] 

We have checked that the Lyapunov exponents in these regions are quite 
different from each other, at least when one integrates the equations of motion 
up to i = 10^. Of course, if an orbit lies on the "boundary" between two of 
these domains, if integrated long enough, it may drift from the inner to the 
outer chaotic region, where its Lyapunov exponents are expected to change 
accordingly. 

In section we will study in more detail the relative location of chaotic 
domains in phase space and argue that these will "overlap" when their respective 
Lyapunov spectra begin to converge as a function of increasing energy for fixed 
TV. 



2.3 Comparison with results in the hterature 

It was shown very recently in that the linear modes of the FPU /3 - model can 
be continued as SPOs of the corresponding nonlinear lattice, i.e. as exact, time 
periodic solutions, having nearly the same spatial configuration and frequency 
as in the linear case. These new solutions arc characterized by exponential 
localization in the q - space of the normal modes Qj{t),j = 1,...,2N and 
preserve their stability for small enough f3 > 0. In fact, the energy threshold for 
the destabilization of the q = 3 solution found by 0| coincides with the "weak" 
chaos threshold determined by de Luca and Lichtcnbcrg in 0. 

In this section, we show that the energy threshold found in Q and ^^so 
appears to coincide with the instability threshold of the SP02 mode. This is 
somewhat surprising since, in all studies of the breakdown of FPU recurrences 
so far, the wave number q of the periodic solutions responsible for the transition 
to "weak" chaos is low (g < 4), while our SP02 mode has a considerably higher 
wave number, q ~ 2{N + l)/3. 
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Figure 3: (a) The period - doubling bifurcation of SPOl for = 11 at the 
energy E = 2.2. (b) The Lyapunov spectrum of two orbits, one starting very 
close to SPOl and another a little further away, for the same N and E as in 
panel (a), (c) The complex instability of SP02 for = 11 at the energy E = 0.5 
after its first destabilization. (d) The Lyapunov spectrum of two similar orbits, 
one in the immediate vicinity and another more distant from SP02 for the same 
N and E as in panel (c). In this figure /3 = 1. 
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Figure 4: The "figure eight" chaotic region for initial conditions in the immediate 
vicinity of SPOl (~ 10^^), a vague resemblance to "figure eight" for initial 
conditions a little further away (~ 10^^) and a large scale chaotic region in the 
energy surface for initial conditions more distant (~ 1) for N = 5 particles, 
when it is unstable, on the Poincare surface of section {xi,xi) computed at 
times when a;3 = 0. In this picture we integrated our orbits up to tn ~ 10^ in 
the energy surface E = 7.4. 



Thus, using different approaches, the authors of Q and |lOl | report an ap- 
proximate formula, valid to order 0(-^), for the destabilization energy of the 
q - breather solution with wave number 9 = 3 given by 

~ —, r. (23) 

In Fig|31 we compare the approximate formula (|23|l (dashed line) with our 
destabilization threshold for SP02 obtained by the monodromy matrix analysis 
of subsection ESI (solid line), for /3 = 0.0315, following [13. We clearly observe 
very good agreement between our numerical results and those of (|23|l . 

The reasons for this agreement are not yet clear to us. Of course, if one 
wanted to look for a chaotic transition via the destabilization of SPOs, it would 
be natural to analyze first SP02 rather than SPOl, since SP02 becomes unsta- 
ble for lower energies. Perhaps the power law cx XjN"^ is important as it is 
the same in formula (|23|l as well as our SP02 stability results. Still, the agree- 
ment between the two curves in Fig. cannot be due only to the coincidence 
of power laws. The proportionality factors in the corresponding formulas must 
also be nearly the same. 



3 Convergence of Lyapunov spectra 

Let us now start from the neighborhood of the SPOl and SP02 modes and 
examine systematically the onset of large scale chaos in the system as the energy 
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Figure 5: The solid eurve corresponds to the energy E2u{N) of the first desta- 
bihzation of the SP02 for /? = 0.0315 obtained by the numerical evaluation of 
the eigenvalues of the monodromy matrix, while the dashed line corresponds to 
the approximate formula (|23(l for the q - breather solution. 

is increased above the instability thresholds Eiu{N) and E2u{N). To do this, 
we need to evaluate the Lyapunov exponents Li, i = 1, . . . , 2N near the SPOs 
(ordered as Li > L2 > ... > L2n)i which measure the rate of exponential 
divergence of nearby orbits in different directions of phase space as time goes 



Note that Fig. |21 reveals that when one calculates Lyapunov exponents 
starting very close (< 10^^^) to a periodic orbit, one finds that they are closely 
related to the eigenvalues of the monodromy matrix of the local dynamics. In 
Fig. OJc), (d), we have plotted these quantities for the SP02 mode at an energy 
where it has undergone a complex bifurcation and has two pairs of eigenvalues 
off the unit circle on the complex plane. 

For comparison, we have also calculated in Fig. Ota), the eigenvalues of the 
monodromy matrix of the SPOl mode at an energy where it is unstable with 
only one pair of eigenvalues off the unit circle on the real negative axis. As we 
see in Fig. Olb), very near this orbit the Lyapunov spectrum is again in close 
agreement with the local results. Of course, in both cases, starting with initial 
conditions a little further away yields the true spectrum, as a smoothly decaying 
curve, which is an invariant of the dynamics in that region (see the dashed lines 



Thus, in the neighborhood of unstable SPOs, one can easily find evidence 
of "small" scale chaos, which is visible at energies where these orbits have just 
destabilized. This, however, is only a local effect, which may have nothing to 
do with the chaotic behavior anywhere else in the system. How could we use 
the dynamics near SPOs, to obtain more global properties of the motion, like 
e.g. the onset of large scale chaos in phase space? 

One way to answer might be to test whether the chaotic regions of the two 
SPOs become "connected" in phase space, above a certain value of the energy. 



to infinity [IS Hi 



in Fig. lab), (d)). 



12 



Evidence that such "merging" of chaotic regions in phase space indeed occurs 
can be provided by their maximal (Li) Lyapunov exponents becoming equal 
and, more specifically, by the convergence of the corresponding Lyapunov spec- 
tra in their vicinity to an exponential function with a characteristic exponent. 

To see this, let us proceed to exhibit in Fig. E^a) the Lyapunov spectra of 
two neighboring orbits of the SPOl and SP02 modes (all orbits in this figure 
start at distances ~ 10~^ from the SPOs in phase space), for iV = 11 degrees 
of freedom and energy values Ei = 1.94 and E2 = 0.155 respectively, where the 
SPOs have just destabilized. As expected, in this case, the maximum Lyapunov 
exponents Li, are very small (« 10"^) and the corresponding Lyapunov spectra 
are quite distinct. 

Turning now to Fig. Elb), we observe that at the energy value E = 2.1, the 
Lyapunov spectra for both SPOs are much closer to each other, even though 
their maximal Lyapunov exponents Li are still different. Furthermore, in Fig. 
Etc), at E = 2.62, we see that the two spectra have nearly converged to the 
same exponentially decreasing function, 

L,{N)(xe-°'TT,i = l,2,...,K{N) (24) 

at least up to K{N) fa while their maximal Lyapunov exponents are vir- 
tually the same. The a exponents of H24|) for the SPOl and SP02 are found 
to be approximately 2.3 and 2.32 respectively. Finally, FiglSJ^d) shows that this 
coincidence of Lyapunov spectra persists at higher energies. 
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Figure 6: (a) Lyapunov spectra of neighboring orbits of SPOl and SP02 respec- 
tively for = 11 at energies E = 1.94 and E = 0.155, where they respectively 
have just destabilized, (b) Same as in panel (a) at energy E = 2.1 for both 
SPOs, where the spectra are still distinct, (c) Convergence of the Lyapunov 
spectra of neighboring orbits of the two SPOs at energy E = 2.62 where both 
of them are unstable, (d) Coincidence of Lyapunov spectra continues at higher 
energy E = 4. All initial distances between nearby trajectories are ~ 10^^. 
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We regard this coincidence as an indication that the chaotic regions of the 
two SPOs have "merged" in phase space in the sense that they "communicate" , 
as orbits starting initially in the vicinity of one SPO may now visit the chaotic 
region of the other. This is strong evidence of the existence of large scale chaos 
in the FPU lattice, at least over the part of phase space travelled by the SPOl 
and SP02 orbits, during their time evolution. 

4 Conclusions 

In this paper we investigated the connection between local and global dynamics 
in the FPU (3 - model from the point of view of stability of its SPOs. Initially, 
we showed that a relatively high q mode of the linear lattice, with one particle 
fixed every two oppositely moving ones, called SP02, is stable for low energies 
until it undergoes complex instability. In parallel, we also studied the properties 
of another mode called SPOl, which keeps every two particles fixed, with the 
ones in between executing exactly opposite oscillations. 

Our first result concerning these orbits is that the energy threshold of the first 
destabilization of SP02 goes to zero faster than that of SPOl. Additionally, we 
discovered that, as a function of iV, the SP02 destabilization threshold coincides 
with the one found by other researchers for the transition to "weak" chaos and 
the destabilization of the 9 = 3 mode. This implies that if chaos is to spread as 
a result of the destabilization of SPOs in the FPU lattice, one might as well look 
closely at the properties of the SP02, as that becomes unstable much earlier 
than SPOl, as A'' increases arbitrarily. 

In order to examine their local dynamics in more detail, we raised the energy 
above the destabihzation of SPOl and SP02 and calculated the Lyapunov spec- 
tra in their neighborhood. We thus found that, as E increases, the Lyapunov 
spectra in the neighborhood of these SPOs appear to converge, at some rela- 
tively low energy value, to the same functional form, implying that their chaotic 
regions have "merged" and large scale chaos has spread in the FPU lattice. 

We, therefore, argue that by studying local dynamics near some of the sim- 
plest periodic orbits which destabilize at low energies, one can gain a better 
view of the "global" dynamics of nonlinear lattices. Furthermore, by computing 
and comparing Lyapunov spectra in the vicinity of such orbits, it is possible 
to obtain valuable insight into the conditions for (or obstructions to) full scale 
chaos and ergodicity, so that we may be able to define probability distributions 
and compute thermodynamic properties of the lattice, as the energy and the 
number of degrees of freedom increase indefinitely, while the energy density is 
kept fixed. 
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